September 24, 2009 16:17 WSPC/INSTRUCTION FILE everitt'2008-09-02 



Non-linear dynamics, entanglement and the quantum-classical crossover of two 

coupled SQUID rings 



M. J. Everitt 

Centre for Theoretical Physics, The British University in Egypt, El Sherouk City, Postal No. 11837, P.O. Box 43, 

Egypt- 
Department of Physics, Loughborough University, Loughborough, Leics LE11 3TU, United Kingdom 

m.j. everitt® physics, org 

Received Day Month Year 
Revised Day Month Year 

We explore the quantum-classical crossover of two coupled, identical, superconducting quantum in- 
terference device (SQUID) rings. We note that the motivation for this work is based on a study of a 
similar system comprising two coupled Duffing oscillators. In that work we showed that the entangle- 
ment characteristics of chaotic and periodic (entrained) solutions differed significantly and that in the 
classical limit entanglement was preserved only in the chaotic-like solutions. However, Duffing oscil- 
lators are a highly idealised toy model. Motivated by a wish to explore more experimentally realisable 
systems we now extend our work to an analysis of two coupled SQUID rings. We observe some differ- 
ences in behaviour between the system that is based on SQUID rings rather than on Duffing oscillators. 
However, we show that the two systems share a common feature. That is, even when the SQUID ring's 
trajectories appear to follow (semi) classical orbits entanglement persists. 
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1. Introduction 

In this paper we report the results of an initial study into the entanglement properties asso- 
ciated with the quantum classical crossover of two coupled SQUID rings. This work can be 
considered as being part of a much larger study spanning many years (see, for example, [UEl 
|3l|4l|5l|6ll7]l8l|^[l0l[TTl[I2[l3][l3[l5l[IS[I3- The overarching theme being the quantum- 
classical crossover of systems the exhibit non-linear dynamics in their classical motion. 
Our interest arises from the fact the Schrodinger equation is linear. Hence, when applying 
the correspondence principle, how can one recover classical like trajectories of classically 
non-linear systems in terms of expectation values of quantum mechanical observables? A 
general solution appears to take the form of introducing environmental degrees of freedom 
together with a suitable expression of the correspondence principle. 

This particular study was motivated by a previous investigatiorU^of two coupled duff- 
ing oscillators. In that work we wishe d to extend the analysis of the quantum classical 
crossover of a sing le Duffing oscillatoiQEESEEMini 

to more complex systems. Specifi- 
cally we wished to study the the recovery of classical like orbits from a quantum description 
of two coupled identical Duffing oscillators. Here our interest lay in understanding how 
the purely quantum mechanical phenomena of entanglement would change as the coupled 
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system approached the classical limit. We showed "that the entanglement characteristics 
of two classical states (chaotic and periodic solutions) differ significantly in the classi- 
cal limit. In particular, we show[ed] that significant levels of entanglement are preserved 
only in the chaotic-like solutions Although the Duffing oscillator has been studied for 
many years with great interest it is, nevertheless, a toy model. Hence, recently, our fo- 
cus has moved to include the study of SQUID rings. These are devices that can manifest 
similar properties as the Duffing oscillator (such as a double well potential and classical, 
dissipative, chaos) but are also readily fabricated. Hence, in this paper we investigate the 
entanglement characteristics for the quantum-classical crossover of a system comprising of 
two identical coupled SQUID rings. 

The correspondence principle in quantum mechanics is usually expressed in the form: 
"For those quantum systems with a classical analogue, as Planck's constant becomes van- 
ishing ly s mall the expectation values of observables behave like their classical counter- 
parts'^^. Applying this expression of the correspondence principle together with intro- 
duction of environmental degrees of freedom via unravelling the master equation is the 
standard method of recovering trajectories of classically non-linear systems. Indeed, we 
used this procedure to recover classical like trajectories for our work on two coupled Duff- 
ing oscillators. However, in[T7]we found that this formulation of the correspondence prin- 
ciple cannot be applied to SQUID rings. Instead we proposed a more pragmatic version as 
follows: Consider h fixed (it is) and scale the Hamiltonian so that when compared with the 
minimum area H/2 in phase space: 

(a) the relative motion of the expectation values of the observable become large and 

(b) the state vector is localised. 

Then, under these circumstances, expectation values of operators will behave like their 
classical counterparts^^]. 

In this work the the evolution of the state vector \ip) of the open system comprising of 
two coupled SQUID rings subject to Ohmic damping will be modelled using a quantum 
state diffusion unravelling of the master equation. This model takes the form of an Ito 
increment equation given bjD^El 

|#) = - l -H S ys\ij)dt 



E 
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where the first term on the right hand side of this equation is the Schrodinger evolution 
of the system with Hamiltonian H sys . The effect of the environment is modelled by the 
addition of a second (drift) and third (fluctuation) term that describe the decohering effects 
of the environment on the systems state vector. In this work we consider only one Lindblad 
operator per subspace Lj = ^/2(dj, where a,j is annihilation operator. The time increment 
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Fig. 1 . The dynamics of the expectation values of position and momentum as a function of time (normalised to 
drive periods) for j3 = 0.01, = 0.1, fi = 0.25 and /3 = 1.0. The dynamics have been taken over the same 
duration and displayed so that, with the exception (3 = 1.0, they approximately align at the time (r e ) at which 
entrainment occurs. (NB: Figure and caption reproduced from 14) 



is dt and the d£ are complex Weiner increments satisfying d£ 2 = d£ = and d£d£* — 
din where the over-bar denotes the average over infinitely many stochastic processes. 

2. Review of the entanglement properties of two coupled Duffing oscillator 

In order to set the scene for our new results we now summarise some of the text and results 
from our earlier work on entanglement dynamics in coupled Duffing oscillators^. In that 
work we extended the analysis of previous work on such system a^^l and considered two 
identical, coupled Duffing oscillators. The Hamiltonian for each oscillator was given by 

tt 1 2 i ft 2 4 1 2 , 9i r,\ . . . 

Hi = -Pi + —q, - -q t +— cos (t) qi + —{q l p % + p l q l ) (2) 

where qi and pi are the position and momentum operators for each oscillator. As usual for 
Ohmic damping, the Lindblad operators were simply Li = \/2Tja i (for i = 1,2), where 
<2j is the annihilation operator. As with previous work in this field m [J4] we set the 
parameters gi = 0.3 and Ti = 0.125 respectively. The Hamiltonian for the coupled system 
then was taken to be: 

H = H l +H 2 + nqiq 2 (3) 

where we set the coupling strength /i = 0.2. In this work the parameter (3 is a scaling 
parameter for h with 0=1 being the fully classical limit. 

Solutions to the equation of motion for the state vector (Q]l are shown in Figure Q] Here 
we show the evolution of the expectation values of position and momentum for the two 



24,2009 16:17 WSPC/INSTRUCTION FILE everitt'2008-09-02 



4 M.J. Everitt 



^ 0.3 

CO 




Fig. 2. Mean entanglement entropy as a function of j3 for the chaotic-like and periodic (entrained) states. Here 
we see that the entanglement entropy for system in the chaotic state does not vanish as /3 approaches the classical 
regime. (NB: Figure and caption reproduced from 14) 

coupled oscillators as a function of time for a four different values of /?. The emergence of 
classical trajectories for small (3 is clearly apparent. Furthermore for j3 = 0.01, 0.1 and 0.25 
chaotic behaviour of the ceases after some time and the motion becomes almost periodic. 

We measured the entanglement between the two oscillators by computing the entropy of 
entanglement for the system^. That is the von Neumann entropv^ f t h e reduced density 
operator pi for one of the oscillators i.e. 

S(pi) = -Tr[pi In pj]. 

We then demonstrated that that when these coupled Duffing oscillators enter the almost 
periodic entrained state the entanglement falls rapidly as the system approaches the clas- 
sical regime. However, when both oscillators undergo chaotic like motion we found that 
significant average entangled remained even when the dynamics of the system appear to 
be classical. We found these results (shown in figure[2]) surprising, as entanglement, which 
does not have a classical counterpart, persisted even for classical looking trajectories. 

3. RSJ model and scaling the dynamics 

As already noted, in[l2]we found that the SQUID ring Hamiltonian was incompatible with 
the correspondence principle phrased in terms of making Planck's constant vanishingly 
small. A solution to this problem was found by scaling the system so that the underlying 
dynamics remain, qualitatively, the same. We can see how this may be achieved for SQUID 
rings by looking at the classical equations of motion that are given by the resistively shunted 
junction (RSJ) model. The RSJ equation of motion for the magnetic flux, $, within a driven 
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SQUID ring iJZD 

d 2 <5> 1# <S>-$ X T . /2tt$\ t , , . 

c it? + r -t + —r- + /c si n^ ) = h sm {ujdt) (4) 

where $> x , C, I c , L and R are, respectively, the external flux bias, capacitance and critical 
current of the weak link, the inductance of the ring and the resistance. The drive amplitude 
and frequency are I d and uj d respectively and <f>o = h/2eis the flux quantum. 

We can then rewrite (0| in the standard, universal oscillator like, form by making the 
following definitions: lo = l/y/LC, t = u> Q t, ip = — <J> x )/$o, fx = ®x/&o, P = 
2ttLI c /& , uj = uj d /uj , ip d = 7 d L/$o an d C = 1/2u)qRC. This yields the following 
equation of motion: 

d 2 ip dip j3 

— T + 2Q— + ip + — sin [2tt (ip + ip x ) = p d sin ur (5) 
dr dr Zir 

In this system of units we then see that we can scale the system Hamiltonian through 
changing either C — > aC or L — > 6L so long as we also make the following changes: 
i? — ► y^b/aR, Id ^ Id/ Vb and — > w<j/ Vafr- 

In this paper we use the following basic circuit parameters C = 1 x 10~ 13 F, i = 
3 x 10~ 10 H, i? = 100ft, f3 = 2, w rf = lu q and 7 d = 0.9 fiA. We note that we have biased 
the ring at the half flux quantum, $ x = 0.5$o> so that the potential approximates a double 
well. We change a so that C varies between 1 x 10 _16 F (quantum limit) and 1 x 10 _9 F 
(classical limit), changing other circuit parameters in line with the above methodology. 



4. Quantum mechanical description of the SQUID Ring 



We now turn to the quantum description of this system. The SQUID ring Hamiltonian for 
each of the identical rings on their own \¥^. 



ft Qi 

Hi ~2C 



($,•-**,(*))' 
2L 



hie 

2e 



■ cos 



2tt$, 



(6) 



where the magnetic flux threading each ring, <3>i, and the total charge across each weak link 
Qi take on the roles of conjugate variables for the system with the imposed commutation 



relation 



= ih and $o = h/2eis the flux quantum. Here $ Xi {t) is the external 
applied magnetic flux and incorporates the drive term for each ring. 

We define dimensionless flux and charge operators in the manner usual for the simple 
harmonic oscillator: Xi = yj Cwq j ffi>i and pi = y/\ /HCuJoQi. 
H[ = Hi/hujQ we find that 



2 



Xi(t)f 



2eu>Q 



cos (ftaij) 



(7) 



where ft 



(4e 2 /S)^/(I7C) 



1/2 
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Fig. 3. Mean entanglement entropy as a function of Capacitance two coupled SQUID rings. Here we see that 
the entanglement entropy for system does not vanish even as it approaches its classical limit. 

We note that this method of introducing Ohmic damping does not bring with it the 
frequency shift that arises through the damping term in the classical dynamics. We re- 
solve this p roblem , as for the Duffing oscillator, by the addition of an extra term to the 

[ToTcl 

Hamiltonian That is, the Hamiltonian 

2 I Zeuo I 

is one we can use to generate trajectories that will be comparable to those predicted by the 
RSJ model. 

H total = V ($+ [ ^~^ )]2 --^co S (0^) + 

c 



(PiXi + Xipi) } + fJ,XiX 2 



where [i = 0.2 



5. Results and comparison with the Duffing oscillator 

We are now in a position to present results for two coupled SQUID rings for compari- 
son with figure [2] Here we have computed the average entanglement as a function of ca- 
pacitance (for a selection of scale parameters a). However we note that the entanglement 
entropies presented here are is the average entanglement over either a long time period 
or many similar trajectories. It is not the entanglement associated with the average den- 
sity operator taken of many experiments. This average entanglement cannot therefore be 
considered usable in a quantum information sense. In figure[3]we show this average entan- 
glement entropy. Here the averaging of each trajectory was determined on a point by point 
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Fig. 4. The dynamics of the expectation values of normalised flux and charge as a function of dimensionless 
time for a range of capacitances. 

basis. A sufficient averaging was used so as to ensure that the results presented here had 
settled to within a percent or so. 

As for the Duffing oscillators, here the mean entanglement does not appear to vanish 
in the classical limit (large capacitance). Another surprising feature in common with the 
Duffing oscillator results is that the average entropy is not maximum at the most quantum 
limit (smallest capacitance). 

From our previous study of the Duffing oscillator as well as that presented in [23] we 
expect this average entanglement of the pure state to be associated with and underlying 
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chaotic or non-linear dynamics. In figure [4] we show sample phase portraits for the ex- 
pectation values of flux and charge as a function of dimensionless time. For very small 
capacitance's we see the motion is dominated by noise fluctuations. As the capacitance in- 
creases we see, as expected, the emergence of non-linear dynamical behaviour. However 
rather surprisingly we see that for large values of capacitance (i.e. the classical limit) this 
non-linear dynamics gives way to a nearly periodic motion. This result for SQUID rings is 
in stark contrast to those for the system comprising of Duffing oscillators. Whilst we find 
these results most interesting we are, here, reporting our preliminary results in the study of 
coupled SQUID rings in the correspondence limit. We will explore these results in more 
detail in later work. 

6. Conclusion 

We have shown that for two coupled SQUID rings the average entanglement persists as the 
rings approach their classical limit. Indeed the trend seems very similar to that associated 
with the chaotic motion of two coupled Duffing oscillators. We note that for the uncoupled 
rings the dynamics remains non-linear for all values of the scale parameter a. However, we 
have seen this is not the case for the coupled system. Unlike the entanglement entropies 
for two coupled Duffing oscillators with the SQUID ring based system we see no sudden 
change in the mean entanglement entropy associated with this change in the underlying dy- 
namics. Specifically, even in nearly periodic orbits the SQUID rings exhibit non-negligible 
entanglement entropies even when the underlying dynamics appears almost classical and 
periodic. As this entanglement is associated with nearly periodic motion it may form the 
basis for extracting usable entanglement from this system. This is something that we will 
explore in a more detailed investigation of this system at a later date. 
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